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ABSTRACT 

We consider luminescence in photo-excited neutral 6*60 using the Su-Schrieffer- 
Heeger model applied to a single C 60 molecule. To calculate the luminescence we use 
a collective coordinate method where our collective coordinate resembles the displace- 
ment of the carbon atoms of the Hg(8) phonon mode and extrapolates between the 
ground state "dimerisation" and the exciton polaron. There is good agreement for 
the existing luminescence peak spacing and fair agreement for the relative intensity 
We predict the existence of further peaks not yet resolved in experiment. 
PACS Numbers : 78.65.Hc, 74.70.Kn, 36.90+f 



Present and permanent address. 



1 



In a recent experiment, Matus, Kuzmany and Sohmen[I[] measured lumines- 
cence from Cqq films and interpreted their results in terms of an exciton polaron. The 
main purpose of this letter is to support the above author's interpretation with some 
simple calculations and to clarify a few points. 

We interpret the luminescence in C 60 within the Su-Schrieffer- Heeger (SSH)@ 
model applied to a single C 60 molecule||. That is, we work with the Hamiltonian 

H = E { -(t-a(X l3 -a)) Y\c\ aCj , CT+ tf. c .] + y (X y -a) 2 }+f/]r n iT n u +l/2 £ m{df 3 jdtf. 

<i,j> cr i j 

(1) 

Here fj is the cartesian coordinate of the jth carbon atom, a is the bare carbon-carbon 
bond length and X i} - is the distance between the ith and jth carbon atoms. Our expe- 
rience with luminescence in conducting polymers Q leads us to two approximations: 
1) We neglect intermolecular hopping (we discuss this assumption later). 2) We set 
U=0 (explicit electron-electron interaction is neglected). When considering vibra- 
tional properties of conducting polymers this is a good zeroth order approximation. 

A number of authors have treated the above Hamiltonian using these ap- 
proximations under the further restriction that the lattice (the 60 carbon atoms) is 
treated classically. For our purposes, the most important result of these studies is the 
formation of the exciton polaron when an electron is promoted from the H u orbital 
to the T\ u orbital. That is, the lattice distorts in the sense that the dimerisation ( 
the difference in the two different bond lengths) is virtually destroyed on a ring cir- 
cling the Ceo molecule. Concurrent with the lattice distortion, two electronic states 
are pulled into the gap. Electrons occupying these states live predominantly on the 
distorted part of the lattice (i.e. the ring). It is important to note that the exciton 
polaron does not break the inversion symmetry of the system. 

Our picture of luminescence in 6*60 is then as follows: By a complex dynam- 
ical process the photoexcited C 60 molecule evolves into the state where one electron 
occupies the lower gap energy level and one electron occupies the upper gap energy 
level. By spontaneous emission the system then decays into the electronic ground 
state and a possibly excited lattice vibrational state. A complication here is that the 
exciton polaron does not break parity and the transition from the highest occupied 
molecular orbital to the lowest unoccupied orbital is dipole forbidden. 

The above picture entails that we must treat the lattice quantum mechanically. 
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We continue to work within the adiabatic approximation. However it is difficult to 
work in the adiabatic approximation without further approximation since the lattice 
has 180 degrees of freedom. We therefore adopt the collective coordinate method Q , 
reducing the problem of 180 degrees of freedom to a single judiciously chosen collective 
coordinate. This method has been used with reasonable success to calculate the 
absorbance and luminescence in nondegenerate conducting polymers [Q. 

The key ingredient for the collective coordinate method is, not surprisingly, a 
good choice for the collective coordinate. In conducting polymers, a good collective 
coordinate has been found to be a one parameter family of lattice configurations 
that extrapolates between the ground state and first excited state classical lattice 
configurations. In the case of Ceo, the ground state of the lattice is dimerised, the 
bonds separating hexagons from hexagons, h-h bonds, have a length l\ and the bonds 
separating pentagons from hexagons, p-h bonds, have a length l?,. Experiment gives 
a value of l\ = 1.40 A and / 2 =1-45A. For the first excited state calculations with 
the SSH model tell us the bond length pattern is largely the same as that of the 
ground state other then on a ring circling the molecule. On this ring, consisting of 
twenty carbon atoms, the difference in bond lengths between the h-h and h-p bonds 
is suppressed. We therefore choose a collective coordinate u so that carbon atoms 
not on the ring are fixed and for f\ on the ring 



where the index % (1 < i < 20) labels the carbon atoms on the ring. Here we treat the 
ring as a chain along the x-axis with periodic boundary conditions. This is a simple 
and reasonable approximation since the quantity entering into the total energy is the 
bond length difference. The parameter a, the renormalised bond length is equal to a 
- average bond length shrinkage, where a is the unrenormalised bond length. (For the 
parameter values we have adopted the h-h shrinkage « .11 A and the p-h shrinkage 
~ .16 A, the average shrinkage being therefore « .14 A). In this choice of collective 
coordinate, u=0 approximates the distortion of the exciton polaron and u= (p-h bond 
shrinkage - h-h bond shrinkage)/4 approximates the ground state dimerisation. We 
emphasize that our collective coordinate is not an unreasonable approximation to 
the Hg(8) phonons|J obtained from microscopic calculations |J. In particular, these 
calculations show that the twenty carbon atoms on the ring are displaced almost 
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parallel to the x-axis. 

We proceed to examine the consequences of our collective coordinate. The 
lattice kinetic energy (eq. 1) then reads in terms of the collective coordinate u 

20 20 
1/2J2 m(dfl/dt) 2 = 1/2 J2 m(du/dt) 2 = l/2M(du/dt) 2 . (3) 

1=1 i=l 

where M=20m is the mass of 20 carbon atoms. With this kinetic energy it is easy to 
write down the collective coordinate Schrodinger Equation 

- ^dVAft* + V (u)ip = Ei> (4) 

V(w), the adiabatic potential energy is computed as the total energy of the SSH 
Hamiltonian for a fixed value of u. In figure 1 we have plotted the adiabatic potential 
energy for the electronic ground state and first excited state. Of course, to obtain such 
a curve we have used parameter values for t,a and K. For polyacetylene, we find t=1.35 
eV, a=7.0 eV/A and K=53 eV/A 2 work rather well to reproduce the experimental 
optical properties and dimerisation. If we use these parameter values for Cqo we obtain 
a slightly too large bond length difference (ie. .06 A ) and a somewhat too large optical 
gap («2.2 eV). We have consequently adjusted a to give the proper dimerisation by 
decreasing a to 6.3 eV/A. Such a value of a reduces the naive gap (see below) to 1.96 
eV in better agreement with the experimental optical gap of 1.9eV Q. By solving the 
collective coordinate Schrodinger equation in the lower adiabatic potential we obtain 
a series of discrete levels separated by about .21 eV. The energy differences between 
these levels should correspond to the energy difference between vibronic peaks in 
luminescence. There is consequently reasonable agreement between our calculation 
and experiment since Matus, Kuzmany and Sohmen[|l|] report prominent peaks at 1.70 
and 1.52 eV, that is an energy difference of .18 eV. According to our calculation, there 
should be additional equally spaced peaks. We attribute the absence of such peaks 
(which we calculate to have smaller intensity, see below) to experimental uncertainties, 
material problems etc.. An earlier measurement of luminescence in Cqq films by Reber 
et al. [0 , in fact, seems to resolve 3 peaks differing in energy by «.17 eV. Preliminary 
results of Iwasa et al.|| also indicate the presence of additional peaks. Our value of 
energy level differences of .21 seems relatively insensitive to parameter choices, for 
example if we let a=7.0 eV/A we get a level spacing of .2 eV. 
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A consequence of our theory is that the energy difference between peaks in 
luminescence should depend on the isotope of carbon present in the C 60 molecule. 
Since the adiabatic potential is very close to harmonic (for the lower curve) the 
energy spacing depends on the mass of the carbon atom m like m~ 1//2 . It therefore 
may be interesting to do experiments on luminescence in Cqq films made using C 13 ||. 

We next turn to the intensity of the luminescence. The intensity of the lumi- 
nescence is proportional to (with an energy independent constant) 

/oo 
i;*(s)i> f (s)Q(s)ds | 2 . (5) 
-oo 

where ipi, ipf are the initial and final vibrational wavefunctions and Q(s) is the 
electronic matrix element for the lattice configuration with fixed collective coor- 
dinate s. The difficulty here is that in the dipole approximation for an isolated 
Ceo molecule (and our collective coordinate) Q(s) is zero. We expect that oxygen- 
impurities, coexisting CVo and/or solid state physics effects (other C 60 molecules) will 
make Q(s) nonzero even in the dipole approximation. Such effects are not simple 
to estimate, fortunately Q(s) is probably only weakly s dependent. Therefore, we 
can treat Q(s) as a constant and pull it out of the integral. In calculating relative 
intensities Q then doesn't enter, we need only consider the quantity 

/oo 
rM^ f {a)d8 | 2 . (6) 
-oo 

A straightforward calculation yields figure 2. In figure 2, the solid circles are our 
calculation, while the solid curve is the experiment of Matus, Kuzmany and Sohmen|]l[ 
and the dashed curve is the experiment of Reger et al. [0. Our calculation is by 
no means in perfect agreement with experiment, it does however seem to be not 
unreasonable especially for such a simple theory (Actually our calculated relative 
intensities agree remarkably well with the experiment of Reger et al). 

Lastly, we consider absorption. Our model for absorption is that the lowest 
vibrational wavefunction in the first adiabatic potential makes a transition to various 
vibrational wavefunctions in the second adiabatic potential. We have plotted the 
absorption in figure 3. The most intense absorption occurs at 1.9 eV. This is in 
apparent agreement with the onset of absorption reported in ref.[l]. 

We have considered luminescence within a simple model. Our model agrees 
well with the existing experimental energy differences and there is fair agreement with 



5 



the relative intensity. We predict that more extensive experiments will see more peaks 
in luminescence and an isotopic shift for the luminescence peaks if C 12 is replaced by 
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FIGURE CAPTIONS 



Figure 1. Adiabatic potential energy vs collective coordinate u. For our parameter 
values u=.0125A corresponds to a dimerised lattice and u=0 corresponds to the ex- 
citon polaron. 

Figure 2. Relative intensity of luminescence vs energy. The large dots are our calcu- 
lation, the dashed curve is an envelop of the experimental data from Reber et al.J/J 
and the solid curve is from the experiment of Matus, Kuzmany, and Sohmen [[!]]. 

Figure 3. Relative intensity of absorption vs energy. The large dots are our calcula- 
tion. 

Note: Figures will be sent by the conventional mail. Please send your request to 
harigaya@etl.go.jp. 
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